Disorder and the elusive superfluid phase of para-hydrogen 
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The possibility that disorder may stabilize a superfluid phase of para-hydrogen in two dimensions 
is investigated theoretically by means of Quantum Monte Carlo simulations. We model disorder 
using a random distribution of scatterers, and study the thermodynamic behavior of the system as 
a function of the scatterer density. Disorder gives rise to equilibrium glassy phases with no superfluid 
properties. Indeed, the propensity for quantum exchanges of hydrogen molecules is reduced even 
with respect to what is observed if the scatterers are arranged as a regular crystal, a physical setting 
that also yields no superfluidity. 

PACS numbers: 02.70.Ss,67.40.Db,67.70.+n,68.43.-h. 
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I. INTRODUCTION 

The inabiUty of observing experimentally a superfluid 
phase of bulk para- hydrogen (p-H2),^ is due to the crys- 
tallization of the system at a temperature T^ ~ 14 K, 
significantly higher than that (T « 6 K) at which phe- 
nomena such as Bose-Einstein condensation (EEC) and 
presumably superfluidity (SF), may be expected to occur 
in the liquid phase. Crystallization takes place because of 
the depth of the attractive well of the potential between 
two hydrogen molecules, approximately three times that 
between two helium atoms. 

The macroscopic property of SF is theoretically asso- 
ciated with the onset of long exchange cycles, involving 
macroscopic numbers of identical particles. As quantum 
exchanges are greatly suppressed in the crystal phase, 
due to particle localization, it is widely accepted that the 
observation of SF of P-H2 hinges on the stabilization of a 
low-temperature liquid phase (although solid order may 
not be entirely incompatible with superflow.i^i^i^i^ii) 
Several attempts have been made-i^ i^°i^^ to supercool 
bulk liquid P-H2 , but the experimental hurdles have thus 
far made it unfeasible to reach a temperature in the range 
where SF might be observed. 

Other experimental avenues have been explored, aimed 
at suppressing crystallization. They include confinement 
and reduction of dimensionality (which reduces the coor- 
dination number). For example, the phase diagram and 
structure of monolayer P-H2 films adsorbed on graphite 
have been studied by various techniques J^ii^iii^ii^ 

One of the most remarkable aspects^^ is that the melt- 
ing temperature r„j of a solid P-H2 monolayer can be sig- 
nificantly less than that of bulk /^H2. Yet, there appears 
to be a limit to what can be achieved by reduction of 
dimensionality alone. A theoretical study^^ of the phase 
diagram of P-H2 in purely two dimensions (2D) has shown 
that the equilibrium phase of the system at low T is a 
triangular crystal, with no evidence of even a metastable 
liquid phase at lower density. Indeed, the system remains 
solid all the way to the spinodal density. The melting 



temperature Ti„ ~ of the equilibrium crystal is 6.8 K, 
i.e., approximately half that of bulk P-H2 but still signifi- 
cantly higher than the temperature at which the system, 
if it remained a liquid, would turn superfluid, estimated 
at '^ 2 K in 2D. Later studie o"'^'^'^^ of P-H2 adsorbed onto 
weakly attractive substrates, predicted only a very slight 
reduction of Tm with respect to the purely 2D case, due to 
an enhancement of the kinetic energy arising from zero- 
point motion in the direction perpendicular to the sub- 
strate (for example, in the case of lithium, arguably one 
of the weakest substrates existing, a reduction of merely 
0.3 K was observedii). This suggests that mere reduc- 
tion of dimensionality is not sufficient to stabilize a fiuid 
phase of P-H2 down to low enough a temperature for su- 
perfluidity to be observed. 

In 1997, Gordillo and Ceperlej*i^ (GC) proposed that 
a significant reduction of the equilibrium density of 2D 
P-H2, and the ensuing stabilization of a low tempera- 
ture liquid phase, could be achieved by intercalating a 
P-H2 fiuid within a crystal of impurities, incommensu- 
rate with the equilibrium crystal structure of pure 2D 
P-II2. Path Integral Monte Carlo (PIMC) simulations of 
such a model system, yielded some evidence of a possible 
superfiuid transition at T ~ 1 K. However, the size of 
the simulated system was very small (4 impurities and of 
order 10 P-H2 molecules). A subsequent PIMC study^^ 
comprising up to ten times as many particles, strongly 
suggests that a) the equilibrium phase of p-^2 in such 
a setting is not liquid, but instead exhibits long range 
crystalline order, and b) the observation of superfluidity 
in Ref. [19| is merely a finite-size effect (this latter finding 
is confirmed here). 

In this work, we explore theoretically the scenario of 
a 2D P-II2 fluid in a disordered environment. Our aim 
is again to determine whether a non-crystalline equilib- 
rium phase, displaying long-range superfiuid coherence, 
may arise as a result of the frustration of solid long-range 
order induced by the random potential. That disorder 
should promote SF is certainly not obvious, and may 
even be counter-intuitive. In fact, disorder is known to 



promote localization, with the ensuing disappearance of 
SF in a system of hard core bosons^^i However, recent 
numerical work has yielded evidence of a possible "super- 
fluid glassy" phase of condensed helium, characterized by 
simultaneous broken translational invariance (and there- 
fore nonzero shear modulus) and SF, with no diagonal 
long-range order. ^^ It is conceivable that a similar phase 
of P-H2 could arise in disorder. 

We investigate the above scenario by means of Quan- 
tum Monte Carlo simulations based on the continuous- 
space Worm Algorithm (WA). Specifically, wc simulate 
a fluid of P-H2 molecules moving in 2D in the presence 
of a static potential generated by a random distribution 
of identical impurities, averaging the results over several 
independent realizations of the external potential. We 
have assumed a specific simple model potential, so as 
to describe the interaction between each scatterer and 
the P-H2 molecules, and varied the overall strength of 
the disordering potential by changing the density of im- 
purities. The main result of this study is a null one, 
i.e., even though we observe exchanges of small groups of 
molecules, they have a local character, i.e., long permuta- 
tion cycles spanning the entire system do not occur; as a 
result, a finite superfluid signal never materializes in the 
low temperature and thermodynamic limits. We observe 
a glassy phase, with broken translational invariance, but 
we also clearly observe that P-H2 molecules locally recre- 
ate the triangular lattice structure associated with bulk 
2D p-H2,~ interrupted by the underlying impurity ma- 
trix. 

Based on this finding, as well as on the null results 
yielded by various simulation works carried out over the 
past few years, we propose that a superfluid phase of 2D 
p-li2 will not be observed, in the absence of some indi- 
rect physical mechanism capable of renormalizing (weak- 
ening) the effective interaction among molecules. 

The remainder of this manuscript is organized as fol- 
lows: Section |TT] offers a description of the model used 
for our systems of interest, including a discussion of the 
potentials and the justifications for the main underlying 
assumptions. Section IIIII involves a brief discussion of 
the computational techniques employed, in addition to 
details of calibration and optimization. The results are 
presented in Sections IIVI FVl and IVIl finally. Section IVIII 
is a summary of the findings and our concluding remarks. 



(M) impurities and hydrogen molecules are regarded as 
point particles, with impurities fixed in space at posi- 
tions Rfc, k — 1, ..., M. The model quantum many-body 
Hamiltonian is therefore as follows: 
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Here, m is the mass of a hydrogen molecule, {rj} (with 
j—l,2,...,N) are the positions of the hydrogen molecules, 
and rij = jr^ — rj\. V is the potential describing the 
interaction between any two hydrogen molecules, and U 
represents the interaction of a hydrogen molecule with an 
impurity. All pair potentials are assumed to depend only 
on relative distances. Most of the results presented here 
apply to the case where the positions of the impurities 
are random, but obviously the above microscopic model 
can describe regular, crystalline impurity arrays (indeed, 
the very same model was adopted for that purpose in 
Refs. HiandHi). 

The interaction V is described by the Silvera-Goldman 
potentialr^ which provides an accurate description of en- 
ergetic and structural properties of condensed P-H2 ^^^ 
The interaction of a hydrogen molecule with each impu- 
rity scatterer is modeled using a standard 6-12 Lennard- 
Jones (LJ) potential, with parameters e = 9.54 K and 
a — 3.75 A as utilized in Ref. [lOJ. We wish to stress 
that, whether or not such a model potential offers an ac- 
curate description of any realistic interaction is largely 
immaterial for our study, as our aim is to study the gen- 
eral effect of disorder; in fact, in principle just about any 
interaction could be utilized, if disorder is produced by 
randomly placed scatterers of variable density. 

In addition to restricting the dimensionality of the sys- 
tem to d=2, the model ([T]) clearly contains important 
physical simplifications, such as the restriction to addi- 
tive pairwise interactions (to the exclusion of, for exam- 
ple, three-body terms), all taken to be central. However, 
here we are primarily interested in examining the funda- 
mental physics contained in a general many-body prob- 
lem described by a Hamiltonian such as ([1]), rather than 
making quantitative predictions about an actual experi- 
mental system - though it should be noted that an actual 
experimental realization may be feasible i^^'^^ 



II. MODEL 



III. COMPUTATIONAL METHOD 



We consider a system of N identical hydrogen 
molecules (composite bosons) in the para nuclear spin 
state, intercalated within a strictly 2D array of identical 
impurities, the entire system enclosed within a simula- 
tion cell of sides L^ and Ly {A=Lx x Ly), with peri- 
odic boundary conditions in both directions. All of the 



Thermal expectation values for quantum many-body 
systems described by a Hamiltonian such as ([1]) can be 
computed by means of Quantum Monte Carlo (QMC) 
simulations. This technique provides numerical estimates 
that one may regard for all practical purposes as as exact 
(at least or Bose systems, of interest here), as results 



are only affected by a statistical uncertainty that can be 
rendered negligible in most cases. In this study, finite- 
temperature estimates were obtained making use of the 
continuous-space worm algorithm, described in Refs. |28| 
andii. 

The only inputs in the simulations are the Hamiltonian 
and the temperature. As with standard PIMC, the cal- 
culation has a finite time step error, which can be made 
sufficiently small by taking the time step short enough; 
for all of the system sizes explored, it was found that 
convergence of the energy estimate is achieved (within 
statistical uncertainty) using a value of the time step 
Tw = 1/640 K~^. Convergence for properties such as 
the superfluid density is reached using a much longer 
time-step, and this property has been taken advantage 
of, where indicated. Based on comparisons of results ob- 
tained from simulations with different values of the time 
step, we estimate our systematic error on the total energy 
per P-H2 molecule (the observable that is most sensitive 
to time-step error) to be of the order of 0.1 K or less (of 
order 0.6%). 



A. Simulations 

All of the simulation results presented below are ob- 
tained starting from an initial configuration of hydro- 
gen molecules randomly dispersed within the impurity 
array. We have studied the physical behavior of the sys- 
tem at sufficiently low temperature, with the aim of ex- 
trapolating the results to the T—0 limit. Specifically, the 
lowest temperature for which estimates were obtained is 
3.125x10-2 K. 

Our simulations for P-H2 in a disordered environment 
are carried out using two different simulation cells. The 
first is a 30x30 A^ square cell, wherein we arrange 
Af =12, 16, and 20 impurities (corresponding to impurity 
densities 17=0.0133 A-^, 0.01778 A-^, and 0.0222 A'^ 
respectively) in a completely random fashion. The sec- 
ond cell, with dimensions 40x34.64 A^, contains Af=16 
randomly placed impurities (corresponding to 17=0.01155 
A^2). For each impurity density, eight different realiza- 
tions are considered (i.e., different random positions of 
the impurities). 

In a slightly different arrangement, we have considered 
a 40x34.64 A^ cell, in which we first arranged M=16 im- 
purities on a regular triangular lattice (details are given 
below), and then applied to each impurity are random 
displacements Dq^i £(0.0,1.0) A, 150.5,2 e(0.5,2.0) A, 
1^0,2 £(0.0,2.0) A, and Di,2 G(1.0,2.0) A. Eight such 
"nearly-periodic" impurity matrices were realized in each 
case (32 total). For each realization, calculations for a 
range of P-II2 coverages were carried out for selected cov- 
erages. 

For each realization and each set of system parameters, 



energetic and structural properties were studied; in gen- 
eral, these remain unchanged below T=2 K. Superfluid 
properties of the system are investigated through the di- 
rect computation of the superfluid density using the tra- 
ditional winding number estimator;'^'^ we as well obtain 
the one-body density matrix, which is readily accessible 
using the worm algorithm. ^^ 

A detailed discussion of the results of these simulations 
is provided in Sections |V] and IVII However, in order to 
provide a useful baseline for comparison, we shall first 
discuss our results for a regular, crystalline matrix of im- 
purities arranged on a triangular lattice, i.e., the same 
system already investigated in Refs. |l^ and l2Cll. 



IV. RESULTS: P-H2 INTERCALATED WITHIN 
A PERIODIC 2D IMPURITY MATRIX 



We have obtained results for the same three system 
sizes of Ref. 20, i.e., the simulation cell contains triangu- 
lar lattices consisting of 2x2 (cell area is ^ = 20x17.32 
A2), 4x4 {A = 40x34.64 A^), and 5x6 {A = 50x51.96 
A^) impurity atoms with 10.0 A nearest neighbor spac- 
ing. In order to assess the robustness of our predictions in 
the thermodynamic limit, we have also performed simu- 
lations for a 10x12 (100x103.92 A^) lattice of impurities, 
at T=4 K and r=0.25 K (the latter using a 4x larger 
time step). We use these results largely to infer/confirm 
structural properties. 

In order to obtain the ground state equation of state 
of the system, we have computed the energy per particle 
e as a function of coverage 6, extrapolating the results to 
the T — > limit. In all energy calculations, corrections 
from the contribution to the potential energy arising from 
particles (and their periodic images) beyond the interac- 
tion cutoff radii were incorporated. Our results (Fig. [1]) 
are identical, within statistical uncertainties, with those 
of Ref. [20. In particular, no significant dependence of the 
energy on temperature is observed for T below 2 K, and 
even the dependence on the size of the system (for all but 
the smallest one considered) is within, or very close to, 
statistical error. Our results are also in qualitative agree- 
ment with those of Ref. [ijl, although estimates offered 
therein are shifted upwardly by as much as a few K (the 
reasons for this discrepancy remain unclear, as merely 
finite size corrections can not account for the difference). 

The equilibrium 2D density 0^ is 0.0381 A-^ which, 
as we shall show below, corresponds to a commensurate 
superstructure formed by the p-il2 molecules. Further- 
more, a stability analysis of the data for e{9) (shown in 
the inset of Fig. [T]) indicates that doping above 9e will 
not result in a homogeneous phase, but rather in the co- 
existence of two (commensurate) phases. 

Figured] shows a snapshot of the P-II2 world lines for a 
system with equilibrium P-II2 coverage 0=0.0381 A^^ at 
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FIG. 1: (Color online) Energy per p-H2 molecule for a regu- 
lar crystalline arrangement of impurities, at the equilibrium 
coverage ^e=0.0381 Kt^ as a function of the temperature, for 
different system sizes. Different symbols show estimates for 
lattices of 4 (circles), 16 (squares) and 30 (diamonds) impu- 
rities. Inset shows the energy per molecule computed at T=l 
K as a function of the coverage 6, obtained for a system of 30 
impurities. 



T = 0.25 K on the 10x12 cell (corresponding to 400 p- 
H2 molecules, precisely 3-|--| times the impurity density). 
The positions of all P-H2 molecules at each one of the 
imaginary-time slices are shown as discrete paths, while 
the fixed impurities are shown as solid circles. The spread 
of the imaginary-time path for each P-H2 molecule gives a 
measure of their quantum delocalization (zero-point mo- 
tion). 

One can immediately identify a trihexagonal tiling of 
P-H2 corresponding to a kagome lattice commensurate 
with, and having three times the density of the impurity 
background (i.e. three P-H2 per centroid of every nearest- 
neighbor impurity triplet). Superimposed over this struc- 
ture is an additional 1/3 coverage of ^-112 molecules ar- 
ranged on a triangular lattice; sites of this triangular lat- 
tice correspond to regions of the total structure having 
four P-H2 molecules per centroid of nearest-neighbor im- 
purity triplets (schematically shown in inset of Fig. [2]). 
As mentioned above it is found, using the procedure de- 
scribed in Ref. |T^ that this phase is incompressible. 

Fundamental theoretical arguments;^ as well as sim- 
ulation of hard core bosons on the triangular lattice^ 
(with which the commensurate phase of P-H2 described 
above bears significant similarities and seems a close 
physical realization) strongly suggest that the commen- 
surate phase of P-H2 found here is not supersolid. In 
any case, in order to assess the superfluid response of the 
system, one can compute the superfluid density, ps, as 
a function of temperature. Accurate scaling of numer- 
ical results obtained for systems of sufficiently different 



and large sizes is of paramount importance in the pres- 
ence of an underlying crystal, with the ensuing loss of 
translational invariance, in order to afford a reliable ex- 
trapolation of the physical behavior characterizing the 
thermodynamic limit. The system of interest here offers 
a chief example of the importance of this aspect. 

For example, the estimates of ps for the 4x4 impurity 
lattice, which are in agreement with the results obtained 
for the 2x2 system, yield a finite value of ps, approxi- 
mately equal to 25%. However, upon increasing the sys- 
tem size to the 5x6 cell, this finite superfluid signal goes 
down to a mere 5%, vanishing altogether for the 10x12 
cell; concurrently, the one-body density matrix displays 
a clear exponential decay with distance, when evaluated 
on sufficiently large size systems. It should be noted that 
the change of behavior observed on extending the sys- 
tem size is not attributable to some loss of efficiency of 
the simulation method in sampling long permutations of 
P-H2 molecules. While this would be certainly an issue 
with conventional PIMC, the continuous-space Worm Al- 
gorithm does not suffer from any such system size limi- 
tation, as shown by a number of calculations carried out 
over the past few years, e.g., for extended defects in solid 
helium, including several thousand particlesi^^iM 

Rather, a large size is simply required in order for the 
actual, non-trivial equilibrium structure of the system to 
be possible, as shown in Fig. [2] Our analysis reveals 
that P-H2 forms a commensurate crystal phase with no 
interstitials or defects. Such a phase, masked for smaller 
system sizes (for which, while the obtained structure is 
artificially frustrated, remnants of the crystalline phase 
are noticeable only a posteriori), is insulating, '^^ with ex- 
plicit calculations of the superfluid density lending con- 
firmation. Thus, in agreement with Ref. |20, we con- 
clude that all previously reported superfluid properties 
are finite-size artifacts, and that P-H2 in this geometry 
is necessarily a regular quantum solid in the thermody- 
namic limit. 

This does not mean that exchanges of P-H2 molecules 
do not take place; however, they are largely local in char- 
acter; long cycles spanning the whole system and result- 
ing in a finite superfluid response, are exponentially sup- 
pressed in the thermodynamic limit. 



V. RESULTS: P-H2 INTERCALATED WITHIN A 
RANDOM 2D IMPURITY MATRIX 

We now turn to describing the results obtained for 
completely random placement of scatterers. For each of 
the eight independent realizations of a random distribu- 
tion of scatterers at a given density, we compute the 
ground state equilibrium P-H2 coverage 9c', in general, 
the value of 6c depends on the specific impurity realiza- 
tion, particularly for smaller system sizes. It is found to 
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FIG. 2: (Color online) Snapshot of a typical configuration of P-H2 molecules at low T at 9 = 9e, intercalated within the 
periodic 2D fixed impurity matrix (described in the text). The positions of all P-H2 molecules at each one of the imaginary time 
slices are shown as discrete paths, while impurities are shown as solid circles. All distances are in A. Inset shows a schematic 
representation of the ensuing phase, inferred from a visual inspection of configurations. Large circles represent impurities, 
whereas smaller circles are P-H2 molecules. Red circles are P-II2 molecules forming the 1/3 commensurate triangular phase, 
which can occupy either one of two equivalent lattices (different color shades). 



fluctuate typically within ^^10% from one realization to 
another. This gives an idea of the local density fluctua- 
tions that occur in the thermodynamic limit, as a result 
of the inhomogeneity of the system. 

The extrapolated ground state energy per particle, 
e(iV), for a disordered arrangement of impurities, roughly 
interpolates between that for bulk 2D p-}i2 (see, for in- 
stance, Ref. 13) and the energy for P-H2 in the pres- 
ence of periodic impurity arrays of the same density^SS as 
should be expected. 

Fig. [3] shows a snapshot of the P-H2 world lines for 



for a specific realization of the disorder with an impurity 
density Q = 0.0133 A'^ (the temperature is T = 0.5 K). 
By visual inspection, we can clearly establish that the 
P-II2 molecules attempt to recreate the triangular lattice 
associated with bulk 2D p-H2r^ interrupted by the un- 
derlying impurity matrix. This is consistently observed 
for all realizations of the disorder, at all impurity densi- 
ties. Moreover, the observed little or no overlap among 
quantum-mechanical delocalization "clouds" of the dif- 
ferent molecules are consistent with quantum exchanges 
being suppressed in this system. 



a system at the equilibrium coverage 6e ~ 0.052 A ^, The above findings are confirmed by an examination of 
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FIG. 3: (Color online) Snapshot of a typical configuration of 
P-H2 molecules at T= 0.5 K for 9=0e, intercalated within the 
2D matrix of randomly placed impurities with f2=0.0133 A~^. 
The positions of all P-H2 molecules at each one of the imagi- 
nary time slices are shown as discrete paths, while impurities 
are shown as solid circles. All distances are expressed in A. 




FIG. 4: (Color online) Ground state P-H2 pair correlation 
functions {g{r)) for equilibrium P-H2 coverage. All functions 
have been uniformly displaced upwards by their y-intercept 
to ease comparison. Results a through d refer to randomly 
placed impurities with decreasing density (from f2=0.022 (a) 
to 0.01155 (d) A~^). Also shown for comparison are g{r) 
for bulk 2D P-H2 (e), as well as for the case of a triangular 
impurity lattice (f), with impurity density ^2=0. 01155 A^^. 



the P-H2 pair correlation function, shown in Figure [4] for 
disordered matrices of different impurity densities. Also 
shown for comparison is the pair correlation function for 
bulk 2D P-H2 at the equilibrium coverage (0=0.056 A^^, 
from Ref. [la ), as well as for the case of a regular impu- 
rity matrix of density f7=0. 01155 A~^. The progressive 
loss of structure at long distance, as the concentration of 
impurities is increased, is obvious; in contrast, the short- 
range features of the pair correlation function are essen- 
tially unaffected, largely reproducing those of the pair 
correlation function for bulk 2D P-H2. 

Analysis of the one-body density matrix reveals a rapid 
exponential decay, with no signature of a Kosterlitz- 
Thouless power-law decay. For each realization, we calcu- 
late a superfluid density of precisely zero, with quantum 
exchanges between P-H2 molecules being strongly sup- 
pressed with respect to the systems with periodic impu- 
rity matrices; system-spanning permutations, necessary 
for the system to feature a finite superfluid response, are 
never generated. All evidence points to this system being 
a "Bose glass" , devoid of both solid order and off-diagonal 
quasi-long-range order. Clearly, no enhancement of Bose 
statistics is brought about by disorder. In order to gain 
additional insight into the effect of disorder, we next 
study a physical arrangement of impurities that inter- 
polates between the regular periodic array of Section HVl 
and the disordered configuration of Section [V] 



VI. RESULTS: P-H2 INTERCALATED WITHIN 
A NEARLY-PERIODIC 2D IMPURITY MATRIX 

For the case where impurities are located at positions 
obtained applying random displacements I?i,2 to regular 
triangular lattice sites (see Sec. IIII Ap . we find on aver- 
age a 13% enhancement of the equilibrium P-H2 coverage 
obtained in Section Hvl 

Structurally, one may infer from the pair correlation 
function (f ) in Figure S] that the P-H2 is in a liquid-like 
or glassy state. Shown in Figure [5] is a typical world-line 
snapshot at T=0.5 K. Upon examination of many such 
"snapshots" , we confirm that the (very limited) interpar- 
ticle overlap observed in Figure [H] is the norm. 

An accumulation of statistics indicate an enormous 
suppression of quantum exchanges, almost entirely ab- 
sent. In every disordered impurity realization examined, 
the superfluid fraction ps is found to be exactly zero. 
The structure and lack of superflow exhibited when -D1.2 
random displacements are applied, are also observed for 
£'0.5,2 and Dq2 random displacements, with no signifi- 
cant physical changes observed in these three cases across 
all realizations. Altogether, these degrees of disorder im- 
parted to the system are sufficient to induce strong par- 
ticle localization, and eliminate even the local exchanges 
found for periodic arrays of impurities in smaller systems, 
as is the case for the systems studied in Section FVl 

Upon reducing the magnitude of the impurity displace- 
ments to Dqi a, we find significant variation from one 
realization to the next, with between a 0% and 4% en- 
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FIG. 5: (Color online) Snapshot of a typical configuration of 
P-H2 molecules at T = 0.5 K for 0=9e, intercalated within the 
aperiodic 2D fixed impurity matrix with density f2=0. 01155 
A~^. The positions of all P-H2 molecules at each one of the 
imaginary time slices are shown as discrete paths, while im- 
purities are shown as solid circles. All distances are expressed 
in A. 



hancement of the equilibrium P-H2 coverage with respect 
to the system having a periodic array of impurities (cor- 
responding to at most 2 additional interstitial particles). 

In realizations for which the particle density is un- 
changed, the same insulating phase found in Section IIVI 
is recovered, albeit with a very slight broadening of ^(r). 

In the realizations where additional molecules happen 
to be introduced (for example, in realizations where three 
impurities that are part of one nearest-neighbor triplet 
are displaced radially outward from their centroid by or- 
der 1.0 A), ps is zero for all T, and the broadening of 
g{r) with respect to that of Section HVl is significant. 

These results are adequate to conclude that the intro- 
duction of disorder which is strong enough to introduce 
additional particles into the P-H2 crystal leads to acute 
particle localization and a suppression of exchanges; dis- 
order weak enough to not change the equilibrium cov- 
erage does not cause a perturbation in the necessarily 
insulating P-II2 crystal structure. 



ties. We performed calculations based a simple model, in 
which the impurities are assumed static and point-like, 
and p-H2-impurity interactions are given by a Lennard- 
Jones potential. 

For P-H2 intercalated within a periodic array of scat- 
terers, the explicitly observed crystalline structure of hy- 
drogen allows us to conclude with confidence that the 
results seen for the largest system sizes are indicative of 
the physics in the thermodynamic limit, and not a man- 
ifestation of algorithmic deficiency. The finite superfluid 
response observed in earlier simulations^ of small sys- 
tems (where the crystal structure could not be exactly 
realized) are thus attributable to finite-size effects. 

From our simulations of P-H2 in the presence of dis- 
ordered arrays of impurities, we conclude that the in- 
troduction of such disorder enhances localization of P-H2 
molecules, suppressing quantum exchanges. Though dis- 
order can give rise to a glassy phase, exchange frequen- 
cies are suppressed with respect to the exchange statistics 
observed in the systems with periodic arrays of impuri- 
ties, in agreement with the conventional view of disorder- 
induced (Anderson) localization in strongly-interacting 
Bose systems. 

Our primary conclusion is that P-H2 intercalated 
within a 2D periodic or aperiodic array of impurities is 
not a candidate for the observation of superfluidity. 

The prognosis, it seems, for observing superfluidity of 
P-H2 in 2D films, is not an optimistic one. Though it is 
possible that incommensuration may permit an enhance- 
ment of disorder and even the formation of a fluid-like 
state, as shown in recent studies of P-H2 and 0-D2 ad- 
sorbed on a substrate composed of graphite pre-plated by 
a single atomic layer of krypton, ^^'■^^ it is doubtful that 
this particular flavor of disorder can enhance particle ex- 
change, and indeed here appears to instead suppress it as 
well. Based on this collection of null results, we conclude 
that a superfluid phase of 2D P-H2 is not likely to be 
observed, in the absence of some indirect physical mech- 
anism capable of renormalizing (weakening) the effective 
interaction among molecules. 

Barring such a mechanism, superfluid behavior in p- 
H2 may be observed experimentally in finite size systems 
only, e.g., clusters i^i^ 
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